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ABSTRACT 



The radiated acoustic field from a point-driven, 
fluid- loaded plate was determined several years ago. 
However, current operational problems force a re- 
examination in the light of certain acoustic anomalies. 

The types of waves that propagate in a solid plate 
are reviewed for their relative importance. 
Experimental models were constructed in an anechoic 
acoustic tank in an attempt to model the aforementioned 
anomalies . 

It is found that the effect of a plexiglass plate 
placed at the boundary of a method-of-images geometry is 
that of a phase delay only. Therefore, mode conversion 
from compressional waves to flexural waves is thought to 
be insignificant. 

A second experiment is described in which a plexi- 
glass plate was driven with a small longitudinal 
vibrator. Experimental radiation patterns suggest that 
of a quasi-piston. Plexiglass is an acoustically slow 
media, therefore the coincidence effect is not possible. 
Further, its absorption of acoustic waves make it a poor 
choice for experiments where wave propagation is 
important . 
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I. INTRODUCTION 



A. THESIS INTRODUCTION 

The experiments on which this thesis is based were 
attempts to physically model a sound source mounted on 
the ice in the Arctic using a sheet of plexiglass and 
piezoelectric transducers in an anechoic water tank. 
Specifically, the sound field radiated into the fluid 
was of interest. 

There are many wave types that propagate in a solid 
and along a solid-fluid interface. Each was examined 
for its relative importance, so that attention could be 
concentrated on those that are more important to the 
radiation of sound into the water from a sound source 
mounted on the ice. The Ewing and Press [Ref. 1] method 
of solution covers wave propatation in ice far from the 
source. Stein [Ref. 2] considered radiation from sound 
sources mounted on the ice. 

Two experiments were performed. The first was an 
investigation on the effect of a plate placed at the 
free boundary of the classic method-of-images (Lloyd's 
mirror) geometry. Waves emanating from a source near a 
pressure-release boundary will develop an interference 
pattern with waves reflected from the boundary that can 
be detected by a receiver. It is hypothesized that a 
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plate placed at the boundary will affect the inter- 
ference pattern. 

In the second experiment performed# the source was 
attached to the ice in an attempt to duplicate the 
"acoustic mode#" named by Stein [Ref. 2# p.100]# of the 
Arctic ice sound propagation. The theoretical develop- 
ment of the chapter uses the Mindlin-Timoshenko thick 
plate equations to explain the effects of coincidence 
and fluid loading on the radiated sound field from a 
point-driven# fluid-loaded plate. 

B . BACKGROUND 

The classical two-dimensional theory of flexural 
waves in elastic bars in vacuo# using Lagrange's 
equation# was solved in 1889 by Lord Rayleigh under the 
assumption that the wavelengths were long with respect 
to the thickness of the plate. Thus# the equations are 
invalid for high frequency waves or for sharp 
transients. There are two effects that were neglected: 
rotatory inertia and transverse shear deformation. 
Later# however# Rayleigh introduced into the equations 
the effect of rotatory inertia# but because transverse 
shear deformation was not included# this correction was 
too large. Timoshenko [Ref. 3] included the effect of 
transverse shear deformation in a study of flexural 
waves of bars . 
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Mindlin [Ref. 4] in a study of flexural waves in 
plates included the effects of both rotatory inertia and 
shear deformation in his now classic paper in which are 
developed the so-called Mindlin-Timoshenko thick plate 
equations. Dyer [Ref. 5] # using the Mindl in-Timoshenko 
equations# showed the effects of transverse shear to be 
important even for thin plates. Skudrzyk et al [Ref. 6] 
developed a theory of flexural waves in a driven plate 
similar to that of the theory of reverberation in room 
acoustics . 

In 1966# David Feit [Ref. 7] calculated the acoustic 
pressure radiated into a fluid by an infinite plate 
driven by a point force# using the Mindlin-Timoshenko 
thick plate equations. His solution used the saddle 
point method of steepest descent to evaluate the pres- 
sure integral. Stuart [Ref. 8] explored the solution of 
the Mindlin-Timoshenko equations with then-new insights 
concerning the so-called leaky waves emanating from the 
plate. He asserted that previous discrepancies in the 
saddle point technique solution (a coincidence frequency 
lower than predicted# and small but unaccounted for 
lobes in the far-field pattern below coincidence) were 
due to neglected poles occurring near the saddle point# 
called leaky wave poles# and that these influence the 
far-field pressure. 
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Creighton [Ref. 9] contradicted these findings, 
demonstrating that "... the far-field contributions 
associated with . . . the leaky waves may be altered 
merely by a change in the integration path ..." and 
that the contributions decayed exponentially with depth 
into the water and did not represent physically meaning- 
ful free modes, except for subsonic surface waves. 

As the theory of waves in a plate developed, so did 
the theory of wave propagation in the ice. The theory 
of propagation of elastic waves in a floating ice sheet 
was given by Press and Ewing [Ref. 1], by considering 
the ice as an infinite plate and using classical wave 
theory. Press and Ewing [Ref. 10] then extended their 
results to include the effects of the air. A further 
discussion of their work as a whole is found in Ewing, 
Jardetzky, and Press [Ref. 11]. 

In a recent study Stein [Ref. 2] used the Ewing and 
Press method of solution, but solved the problem with 
the aid of modern computer algorithms . The bulk of the 
study was that of the propagation speeds of the various 
modes. His data taken from Arctic and freshwater lake 
ice experiments corraborate with the predicted transient 
arrival times and propagation losses. 

C. WAVES IN A SOLID 

There are several wave types that may propagate in a 



15 



solid and others that occur at a solid-fluid interface. 
A study of the applicability of each type found in the 
literature was done and is included here to give the 
reader a sense of wave types that are relevant to the 
problem of modeling a sound source mounted on the ice 
generating a sound field in the underlying water by the 
excitation of the ice. 

Propagation of waves in solids can be distinguished 
from that in fluids by the solid's ability to restore 
shear. This leads to two wave equations: 



for the scalar displacement potential which describes 
longitudinal motion, and 



for the vector displacement potential which describes 
shear wave motion. The potentials <|>(x,y,z,t) and 
4>(x # y,z,t) are related to the particle velocity vector 




P 



(1.1a) 




(1.1b) 





( 1 . 2 ) 
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An excellent depiction of some of the different 
types of acoustic waves, along with their speeds, is 
given by Ver and Holmer [Ref. 12, p. 273] and it is from 
this source that Figure 1.1 is taken. 

Longitudinal waves are waves that propagate along 
the direction of compression and dilatation in the 
plate. The thickening and thinning of the plate seen in 
Fig. 1.1a is a consequence of a non-zero Poisson's ratio. 

Shear waves are solutions to the additional wave 
equation that governs waves in solids : 



where 




c, 2 - — as in Fig. 1.1b 
s P 



+ 

= vector potential of particle displacement 
in a solid. 



Two types of shear waves are found in the 
literature: SH and SV waves are the designations given 
them in texts on seismological wave propagation. SV, 
shear waves polarized vertically, denotes waves for 
which the shearing motion is perpendicular to the plane 
of the plate. It is this wave that is depicted in Fig. 
lb. SH, shear waves horizontally polarized, denotes 
waves for which the shearing motion is parallel to the 
plane of the plate. 
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shear wave 
C S = 





torsional wave 

c<j» = 



GK 

Pi 



bending (flexural) wave, 
(thin plte theory: see 
text) 

cb 55 



V (i> 2 Yt2 

12 p ( 1- a z ) 



Figure 1.1. Wave Types in Solids and Their Speeds 

(After Ver and Holmer) 

Y * Young's modulus, which relates the stress to the 
strain produced, 

a = Poissons ratio, which relates lateral expansion to its 
longitudinal shortening due to a compressive stress, 
G * the shear modulus, which relates shearing stress to 
shearing strain produced, 

I =* the moment of inertia, 

K =* a torsional stiffness factor, which relates a twist 
to the shearing strain produced, 
p =* the solid's density, 
t = the solid's thickness, and 
to = the radian frequency of the motion. 
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Torsional waves. Fig. lc, result from a twisting 
force applied to a solid. This wave can be generated in 



bars, where the y and z dimensions are similar. 
Torsional waves are not be easily excited in a plate; 
they will not be further considered in this study. 

Flexural waves, or bending waves, are a more complex 
motion than the previous ones considered, because of the 
complicated stresses set up when a plate is bent. The 
bending of a plate compresses the plate on the inside of 
the bend and expands the plate on the. outside. Due to 
this compression and expansion, a Poisson's ratio effect 
is introduced . 

As with any wave, the energy of a flexural wave is 
composed of two parts, potential energy and kinetic 
energy. The potential energy is made up of two 
deflections: flexural deformation, and shear 
deformation. The kinetic energy also has two 
contributions, linear (or translatory) inertia, and 
rotatory inertia. 

The classical two-dimensional theory of flexural 
motions in thin elastic plates leads to the wave 
equation [Ref. 4, p. 31] 




(1.3) 



where D 



Yh 3 



and here 



12(l-o*) 
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h = plate thickness 
w = the plate particle motion 
This equation relates only the flexural deformation to 
the linear (or translatory) inertia as the two opposing 
forces . The Mindlin-Timoshenko equation relates the 
linear and rotatory inertia as the driving forces to the 
flexural and transverse shear deformation as the restor- 
ing forces. (This was what was referred to earlier as 
"correcting" for rotatory inertia and transverse shear.) 
The Mindlin-Timoshenko equation is 







w + 




(1.4) 
w = 0 



where 



r i _ xf. Y 
G 12 TTI+o 7 ) 

h = plate thickness. 



The 



32 

at 7 



term is the contribution of transverse shear. 



The 



ph 3 3 2 
12 . 3t 



term is the contribution of rotatory 



inertia . 



Note that if these terms are neglected, the classical 
equation (Eq. 1.3) results. 
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The Mindlin-Timoshenko theory is used when analyz- 
ing "thick" plates / i.e., plates whose thickness exceeds 
X s /20 , where A s is the shear wave length. [Ref. 13, p. 
215] 

Of key importance to this study is the effective 
phase speed for flexural motions as a function of 
frequency. Ross, assuming a harmonic disturbance of a 
homogeneous plate in vacuo, gives an equation for the 
flexural wave speed Cfc, as 



(i- r p (|-)) 2 (|V 

*'p p 



+ < r P +1 ><%> <=;> 



(■g-) 2 = o (1.5) 
“p 



where 



Tp = the plate shear parameter 
= 2.65 (1 + | a + | a 2 ) 
flp = the reference angular frequency. 




a) = the angular frequency of the flexural wave 

Cp - the plate longitudinal wave speed, given in 
Fig. 1.1a. 

Equation 1.5 can be solved by applying the quad- 
ratic equation to solve for the variable (■^ 2 ) . It is 

C P 

a difficult function of the other parameters, but can be 
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represented accurately in terms of approximate solutions 
at low, intermediate and high frequencies. 

At low frequencies, the flexural wave speed is 
approximated by 



c bl 



w 2 Yt 2 
12 p ( l-a z ) 



( 1 . 6 ) 



which is that given in Fig. l.ld, and represents the 
flexural wave speed in thin elastic plates, given by the 
classical wave equation. Equation 1.3. This equation is 
valid when 



T ( — ) 



2 



<<1 



In the intermediate frequency range, the flexural 
wave speed is approximated by 



°bi “ c bl 




(1.7) 



In the high frequency range, the flexural wave phase 
velocity predicted by the Mindlin-Timoshenko theory 
tends to the Rayleigh wave velocity [Ref. 13, p. 248], 
which is a value slightly less than the shear wave 
velocity. The expression for the high frequency 
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flexural wave speed is 



c bh * 0.76(1 + | o) c s 



( 1 . 8 ) 



This equation is valid when 





>>1 . 



Heretofore we have discussed waves propagating in a 
solid. Waves can also propagate along the boundary of a 
solid-fluid interface in much the same manner as do sur- 
face waves on the ocean. These waves are Rayleigh waves 
(for a half-solid). Lamb waves (for a plate), Stonely 
waves. Love waves, and Crary waves. 

Rayleigh waves ". . .represent elastic perturbations 
propagating near the free boundary of a solid and decay- 
ing with depth” [Viktorov, Ref. 14, p. 1]. Figure 1.2 
is a depiction of the instantaneous deformation associ- 
ated with a Rayleigh wave. 




Figure 1.2 Depiction of a Rayleigh Wave (after 
Viktorov) 
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Rayleigh waves propagate along the surface of a 
solid at a free boundary, such as an air-ice boundary or 
an air-plexiglass boundary. The wave decays exponen- 
tially in the z-direction into the solid, and if the 
decay constant is not large, the motion of the wave may 
appear at the lower (in this case, water) boundary. 
Furthermore, the Rayleigh wave is supersonic with 
respect to the fluid, therefore this motion would radi- 
ate sound into fluid (see Chapter III for an explanation 
on the importance of relative sound speeds to radiation 
into the fluid.) 

The discussion of Rayleigh waves is included 
because, as mentioned above, in the high frequency limit 
flexural waves degenerate to Rayleigh surface waves for 
a plate in vacuo. However, Stein shows that at low fre- 
quencies there is no coupling of the Rayleigh wave with 
the flexural wave there would be for a plate in a 
vacuo 1 . Moreover, the types of sources being considered 
would not generate Rayleigh waves, hence they will not 
be further considered. 

Lamb waves ". . . refer to elastic perturbations 
propagating in a solid plate (or layer) with free 



^There is a growing debate among the cognoscenti 
concerning the importance of the Rayleigh wave in Arctic 
ice [Ref. 15 ] . 
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boundries, for which displacements occur both in the 
direction of wave propagation and perpendicular to the 
plane of the plate." [Ref* 13# p. 673. They are also 
called the "normal inodes" of a plate# and can be sym- 
metrical (s) or antisymmetrical (a) relative to the 
center plane of the plate. Note that the zero-order 
symmetric and antisymmetric modes are the longitudinal 
and flexural waves, respectively. Figure 1.3 depicts 
the zero-order Lamb waves. These two modes are of key 
importance to this study. 




i 




* 



Figure 1.3 Lamb Wave Propagation (after Viktorov) 

Another wave that can exist at an interface between 

two media is a Stonely wave. According to Ewing# 

Jardet.zky and Press [Ref. 11# p. 1123 they are modified 

Rayleigh waves that exist when the shear wave velocities 

• 

in the two media are nearly equal. Stein [Ref. 2# p. 
69] states that they can also exist along the boundary 
between an elastic medium and a high density# 
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acoustically slow fluid or elastic medium/ and that the 
Stonely wave is the high frequency asymptote of the 
first anti- symmetric mode along the ice-water interface. 

That is / they are Rayleigh waves along the solid-fluid 

■> 

interface. However/ absorption in Artie ice is high at 
the frequencies for which Stonely waves exist [Ref. 2, 
p. 190] and therefore the Stonely wave is not important 
to this study. 

Two other wave types concerning propagation of sound 
in ice are described in the literature. Love waves and 
Crary waves . Both are shear waves . Love waves are SH 
waves that are generated due to discontinuities in the 
ice. Stein [Ref. 2 p. 100] notes that "... unless 
there is a discontinuity .... they will not radiate 
nor can they couple to other modes." Discontinuities in 
the ice are not considered in this study, therefore Love 
waves will not be considered. 

Crary waves are a type of SV wave discovered by 
Crary in seismic experiments on a floating ice island. 
They were generated after an ice "event"/ a small ex- 
plosive charge mounted in the ice. According to Ewing, 
Jardetzky and Press [Ref. 11, p. 299], this wave has the 
following characteristics: 

"1. phase velocity is near the speed of compression- 
al waves in the plate; 
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2. travel time is intermediate between that for P 
( compress ional ) waves and that for SH waves? 

3. principal [amplitude] is [observed] on longitud- 
tudinal horizontal seismograph? amplitudes are 
much smaller on vertical and practically absent 
on transverse horizontal? 

4. frequency is almost constant, increasing very 
slightly with time." 

Stein [Ref. 2, p. 93] again discusses these waves, 
but has found the attenuation to be at least 0.1 dB/m 
in the horizontal for the Arctic. Thus, these waves 
will not be observed at appreciable distances. Further, 
because of their small vertical amplitude and velocity 
less than that of the fluid, (characteristic 3 of above) 
they do not impart acoustic energy to the underlying 
fluid. For these reasons, Crary waves will not be 

further considered. 

So far in this chapter, the types of waves that pro- 
pagate in a solid, namely Arctic ice, have been dis- 
cussed without regard to their origin. That is, they 
have been .assumed to be plane waves, far from their 
source. As these waves propagate, they may be coupled 
into the surrounding medium, creating acoustic 
radiation. 

Additional radiation may be generated when the waves 
are excited by a localized force. This is a result of 
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the flexural motion having a spatial gradient of dis- 
placement in the ice, which plane waves do not. Stein 
[Ref. 2 p. 96 J has named this mode of acoustic radiation 
the "acoustic mode." 

Figure 1.4 is a depiction of the first antisymmetric 
mode, the flexural wave, in a solid. At this point, far 
from the source, it can be seen why the flexural wave 
does not radiate 2 : it travels subsonically with respect 
to the water and there is cancellation of radiation from 
adjacent positive and negative pressure peaks. (In 
actuality, because of geometric spreading, there is in- 
complete cancellation and some radiation does exist.) 
However, When these waves are excited by a localized 
force (Fig. 1.5) there is incomplete cancellation and 
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Figure 1.4 Flexural Wave Propagation 



2 A flexural wave below coincidence will not radiate. 
The coincidence frequency is that frequency above which 
the plate vibrations can couple to the fluid. See 
Chapter 3 for a further discussion. 
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radiation occurs. Note that this acoustic mode is still 
a flexural wave. A further discussion of this mode found 
in Chapter III. 
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Figure 1.5 Flexural Wave Propagation From a Point 
Source (after Stein) 

D . SUMMARY 

In this chapter, the many types of waves found in a 
solid (ice) have been discussed. As was indicated 
earlier, only a few have real importance in the Artie, 
and therefore modeling in the lab should be attempted 
only for these. Stein [Ref. 2 p. 107] has concluded 
that the three modes of propagation of importance are 
the first symmetric and antisymmetric modes, and the 
acoustic mode. All other aforementioned waves are 
neglected because either they are highly attenuated in 
the Arctic at the frequencies of concern or there is a 
lack of a source with which to generate them. 
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II. A LLOYD'S MIRROR EXPERIMENT 



A. INTRODUCTION 

In this chapter, an experiment will be discussed 



duced at an air-water interface will influence the 
acoustic field in the water half-space. More simply. 



the reflection of a spherical wave from an ideal 
boundary develops an interference pattern with the 
original incident wave. It is hypothesized that the 
plate will effect this Lloyd's mirror pattern. 

B. THE METHOD OF IMAGES 

The method of images is a well-known technique for 

/ 

predicting the radiation pattern of an acoustic source 
near a fluid interface. For continuity, a brief de- 
scription of the theory of the acoustic Lloyd's mirror 
will be presented. The following description is based 
on that given by Sanders and Coppens [Ref. 16]. 

If a source of spherical waves is placed on the z- 
axis a distance d below a pressure-release boundary, as 
shown in Fig. 2.1, the acoustic radiation from this 
source in the absence of the boundary can be specified 



that tested the following hypothesis: a plate intro- 



t, 



as 



£ s = - e j ( u)t-kr_ ) 



( 2 . 01 ) 



where 
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As a result of the boundary, there must be an addi- 
tional wave present to describe the reflection: 



£r 



- e j(oit-kr + ) 

r+ 



( 2 . 02 ) 



where r + = ^(h+d ) 2 + R 2 



The total field is, then: 



(2.03) 



£s - £r = £ = Mi- e^ kr - - ^ 



The acoustic field will show strong interference effects 
because the two waves £ s and £ r must travel along 
different paths to reach the receiver. 

The geometry is simplified if the assumption 



r>>dcos0 (2.04) 

is made. Then, r+ and r_ can be replaced by 



r + * r+Ar (2.05) 

r_ » r-Ar 

Ar - dsinQ 

with 2.04, 2.03 becomes 



£(r, 0 , t ) 



A j ( o)t-kr ) rej^A r 

r S i Ar 



e-j kAr 
Ar 
r 
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Figure 2.1 Lloyd's Mirror Geometry 

Thus when a simple harmonic acoustic wave in the the 
water is incident on a plane, perfectly reflecting pres- 
sure release boundary, the reflected wave is a mirror 
image of the incident wave which leads to an interfer- 
ence pattern in the medium below the boundary. 

This interference pattern may be observed by dif- 
ferent methods. The receiver can be placed at a fixed 
horizontal range from the source and moved downward to 
show the nodes and antinodes from the pattern. It might 
also be placed at a fixed slant range from the source 
and swept through an arc about the source. Or, it might 
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be held at a constant depth and moved away from the 
source in horizontal range. Lastly, and want was done 
in the experiment, the source and receiver can be held 
in a fixed position and the frequency of the source 
swept, observing the response of the receiver. The 
interference of the incident and reflected waves produce 
maxima and minima in the frequency response of the 
receiver . 

C. THE EFFECT OF THE PLATE 

The effect of the plate 3 on the incident wave can be 
best seen by reverting to simple ray theory. Consider 
the following problem: an incident wave in water, re- 
presented by ray impinges on a plate of different 

characteristic impedance r , as shown in Fig. 2.2. The 

2 

ray (now £ a ) is refracted as it traverses into the 
plate. It reflects from the top of the plate at z = H 
(now , traverses the plate and then is refracted back 
into the water as p r • Note that while the wave is in 
the plate, it is considered as a compressional (longitu- 
dinal) wave . 



3 An infinite plate forced by an incident sound wave 
at any frequency, even below its coincidence frequency, 
will radiate sound, although a bending wave of the plate 
itself, below coincidence, will not. [Ref. 12, p. 280] 
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Figure 2.2 Ray Traversing Through a Plate 



The four waves may be represented as (considering only 
two dimensions): 

Ei = p i ej( wt " k sin0^x-k ^cos 0^z ) (2.1) 

- » P a ejUt-fc sin0 a x-k cos0 a z) (2.2) 

xia — a 2 2 

Eb = Zb ej( “>t-k^sin0 b x+k^cos0 b z) (2.3) 

a e j(u>t-k sin0 r x+k ^cos 0 r z ) (2.4) 



The boundary conditions are: 

at z = 0, continuity of pressure requires 



£ i (x # o) +£ r (x,o) -£ a (x,o) +£ b (x # o) 



(2.5) 
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Also at z=0, the condition of continuity of normal 
velocity, requires (taking into account direction of 
propagation) 

ui(x,o)cose r - UyCx/Ojcosey. 

( 2 . 6 ) 

= U a (x / o)cos0 a -u t) (x,o)coset ) 



The boundary at z = H is pressure release, therefore 

£ a (x,H)+£ b (x,H)=0 (2.7) 

For each plane wave, we have 

£ = ±r (2.8) 

u 

With Snell's Law, s A n 9 r - S A , ^ - 9 - a # (2.9) 

c 1 c 2 

and since the angle of incidence equals the angle of 
reflection, we have 



0j_ = 0 r = 0 1 and 0 a = 0 b = 0 2 (2.10) 

If each value of u in 2.6 is replaced with the 
appropriate value of p/r , and from 2.10, subscription 
are changed, then: 



■ ^ S91 [£i(x,o)-£ r (x,o) ] = [£ a (x,o)-£ b (x,o) ] (2 * 12) 

r l r 2 



or 



£i ( X , O ) -£ r ( X , O ) 



r i cos 0 2 
r 2 cos©! 



(£ a (x,o)-£ b (x,o) ) 



(2.13) 
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Apply 2.1 through 2.4 to 2.5 , and drop the e^ a)t , 



-jk sine x -jk sine x 

(£i+P r ) e 1 1 = ^£a + £b )e 2 2 (2.14) 



and do the same to 2.13, 



( p p )e -i k l sin9 l x = £i£°s0Z (P a -P b )e-i k 2 sin0 2 x (2.15) 
i — r j- cos e — “ — ° 

2 1 



Apply the boundary condition at z=H (Equation 2 . 1 ), and 



-jk sine x 

dividing through by e 2 2 gives 



p^e" jk 2 cos 9 2 H + p- b e + ^ 2COSQ2H = 0 

l 



(2.16) 



Algebraic manipulation of 2.14 - 2.16 then leads to the 
reflection coefficient R : 



R = 



£r 



- £i 



fl - £l cos9 ? 1 _ 
^ r 2 cos 6 t [ 



(1 + -1 
^ r 



cos e 



cos 6 
2 1 



4 ) - (i - 7 



r ^ + r_^ cose? a e " j2k 2 COS 6 2 H 

2 r 2 cos e i 

r , cos 8,^ - j2k 9 cos 9 9 H 



cos e 
2 1 



2 W ~" u 2 
(2.17) 



where e and 9 are related by Snell's Law (.9), and 
l 2 



cos 92 



-aFF 



C ^- 2 ) 2 sin20 i 



(2.18) 



I 



36 



D. VALIDATION OF EQUATION 2.17 

In the limit as H+O , Equation 2.17 reduces to 
r =-l . This is to be expected, because when H+0 the r 
layer disappears and the incident wave is 2 reflected 
off a purely pressure-release boundary, where the re- 
flected wave is 180° out of phase with the incident, and 
hence R = -1 . 

In the limit as 0 >0 , that is, normal incidence, 

1 

Equation 2.17 reduces to: 



Si = 



ri cosk2H + jr?_sink2H 
r ]_ cosk 2 H + jr 2 sink 2 H 



(2.19) 



Comparing Equation 2.19 with Equation 6.12 of . 
Kinsler et al [Ref. 17] 




( 2 . 20 ) 



where r^ is characteristic impedance of the fluid above 

the place. In the limit r +0 , that is, assuming a 

3 

vacuum above the plate, 2.20 reduces to 2.19. 

Finally, comparison of 2.12 with Brekhovskikh [Ref. 
18] yields additional validation. His coordinate system 
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~ Oil 



9h 



med 









zrid 





3a 



which reduces to [Ref. 18, Equation 5.18] 



R 



V? 2 + V i 2 

1 + V V 
23 



q Z j a 2 d 
Z j a 2 d 

12 e 



(2.24) 



where V 12 



z i~ z 2 
Z +Z 
1 2 




Z?-Z 3 

Z +Z 
2 3 



(2.25) 



Assuming Z =0 (in vacuo), then V = -1 , and .27 
1 12 



reduces to 



R 



(1 

(1 




(1 + 

(1 - 




Z^ zja 2 d 
Z ; 



2 



(2.26) 



changing to previous notation and noting that d = -H , 
equation 2.26 becomes 2.17: 



R 



r 1 cos 6 2 1 M + £_l cos 6 ■? 1 

r 2 cos 6 i [ \ £2 cos 6 j ' 

r 1 cos 6? 1 _ f 1 - £j_ cos 9 ? ^ 

r cos 0 ' ^ r cos 6 ' 

2 1 2 1 



e ~z jk 2 cos 0 2 H 
-z jk 2 cos 0 2 H 



E. VALIDATION OF THE PLANE WAVE ASSUMPTION 

Thus far, the reflection coefficient of an oblique 
incident plane wave off a finite layer has been calcu- 
lated and validated against previously known theory and 
limiting cases. However, one question still needs to be 
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addressed before this equation is put to use. That is, 
the validity of the plane wave assumption made above. 
Equation 2.8, needs to be established. 

The solution for a point source near a boundary is 
addressed in two works. Officer's [Ref. 19] solution is 
in the form of three integrals for each of the three 
regions. However, the integrals are difficult to 
compute. Officer does give "... the condition under 
which the plane wave reflection coefficient may be taken 
as a valid approximation for the reflection of an acci- 
dent spherical wave. This condition is that 




1 " 



(2.27) 



where is the distance from point of reflection to the 
receiver . 

The second work that addresses this problem is 
Brekhovskikh . From Brekhovskikh , the geometric optics 
approximation (in which a spherical wave is relected off 
a layer with the same reflection coefficiient as a plane 
wave) is valid when the distances of the source and the 
receiver from the boundary are large compared with the 
wavelength. [Ref. 18 p. 256] In the experiment run, 
this assumption is valid above f = 50kHz. 
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Lastly, for completeness, the extent of the near 
field should be calculated. From Ziomek [Ref. 20], the 
maximum radial extent r of the near field of the source 
is given by: 



r 




(2.28) 



where 

R is the maximum radial dimension of the transmit 
transducer . 

For the experiment done here, the following values 
were used in equations 2.27 and 2.28: 
c^ = 1481 m/ s 
C 2 = 1992 m/s 
0! = 42.27° 

Rj_ = 0.5 m 

Equation 2.27 requires that f >> 14.8 kHz for the plane 
wave assumption to be valid. From 2.28, for the radia- 
tion from the transmitting LC-10 to be in the far-field 
prior to reaching the plate requires f < 1.01 MHz. 

Therefore for the frequency range of interest (50-100 
kHz) both conditions are satisfied. 

F. APPLICATION OF EQUATION 2.17 

Upon close inspection of Equation 2.17, it can be 
seen that the reflection coefficient is a phasor of unit 
magnitude with a phase delay that is a function of 
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frequency. This is intuitively true because of the con- 
ditions imposed on the problem: the boundary condition 
at the top of the plate, Z = H , is a pressure release, 
perfectly reflecting plane; the medium is assumed to be 
lossless. Since it is perfectly reflecting, there is no 
refracted wave at the upper boundary and there is no 
loss of signal. If the medium is lossless there is no 
loss of signal as the wave travels through the plate. 
Thus, the reflection is of unit magnitude with a phase 
angle, that is R = l*e^ T . 

The phase angle x of R is, from Equation 2.17, 



x = tan 



_ x U+-f^ fffff)sin[2k 2 Hcos6 2 ] 



(1- •£->- 52®iz)_(i+ £l ^i2)cos[2k 2 Hoos9 2 ] 
^ r cose ; ^ r cose ; z z 



(2.29) 

[1 - — cos jsinf 2k 2 Hcos6 2 1 
- tan" 1 r? cos9 ' ^ 

i 1 + f 1 - cos 1 2 ) COS C 2k 2 Hcqs 9 2 ) 

2 1 



In the experiment to be described, the following para- 
meters were used in Equation 2.29: 

r l = Pi c i = C998 ^f] [1481 m/s] = 1.48 x 10 6 - a s - 
x A 1 m J m 

r 2 = P 2 C 2 = [1200 *3] [1992 m/s] = 2.39 x 10 6 

m J m 

H = 3.18 x 10 -3 m 
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From the geometry of the setup. 



0 ! = tan" 1 



l rcvr 



T“a: 



“ tan 1 “35m - 42,3 ° 



d rcvr ' “source 
where d rcvr and d source are the receiver and source 

depths, respectively. A plot of the phase angle caused 
by the plate. Equation 2.29, as a function of frequency 
with the parameters shown in Figure 2.4 above is in 
Figure 2.4. Note the conventional of a negative angle 
producing a phase delay in the signal . 

In the frequency range of interest, f = 0 to 100 
kHz, t decreases from 0° at f ■ 0Hz, to -104* at 



f = 100,000 kHz. 



A 




Frequency (kHz) 



Figure 2.4 Phase Delay Caused By the Plate 
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The phase delay caused by the plate will be mani- 
fested by a shift in the frequency of the nodes and 
antinodes of the interference pattern. It can be shown 
that the frequency shift Af caused by the phase delay x 
is : 

4f =-57IF (2 - 30) 



where 

Ar = difference in path lengths between the direct 
and reflected signals 

x = phase delay caused by the plate, as obtained 
from Equation 2.29, 
c 0 = speed of sound in the fluid. 

G. DESCRIPTION OF EXPERIMENT 

Two small LC-10 transducers were placed in the 
water, one acting as a source, the other as a receiver, 
and were connected to a Hewlett Packard 4162A frequency 
analyzer as shown in Figure 2.5. A periodic chirp sig- 
nal of 0-100 kHz was applied to the source transducer, 
and the received response was measured and plotted. The 

i 

LC-10 transducers were assumed to have a flat response 
throughout the frequency range. While the response 
curves of the LC-10 show the response to be fairly flat 
up to 100-kHz, the assumption is really not tested, 
since the response of the transducers themselves 
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is not being tested, rather the difference in response 
with and without the pressure of the plate. 

The response of the receiving transducer showed the 
maxima and minima in the pressure at its location due to 
the boundary. The response was measured for the water- 
air boundary, and without changing the location of the 
two transducers, the plate was introduced and the 
response was remeasured. 

H. EXPERIMENTAL RESULTS 

The results of the experiment can be seen from 
Figures 2.6 and 2.7. Figure 2.6 is a plot of the fre- 
quency response of the receiving hydrophone with fixed 
geometry with no plate present; the classic Lloyd's 
mirror experiment. Note the frequency pattern of the 
nodes and antinodes. When the plate was placed in posi- 
tion without moving source or receiver, the positions of 
these nodes and antinodes shift to the right, higher in 
frequency, as can be seen in Figure 2.7. This was pre- 
dicted as the manifestation of the phase angle caused by 
the plate . The change in the frequency of the nodes 
increases linearly with frequency, corresponding to the 
increased phase angle. 
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Figure 2.5. Lloyd's Mirror Experiment set-up 
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Table 2.1 shows a comparison of the predicted nodal 
shift to the nodal shift obtained experimentally. t was 
obtained from Figure 2.4. The predicted frequency shift 
was obtained from Equation 2.30. 



TABLE 2 . 1 


f node 


T 


^predicted 


exp . 


37 kHz 


0 

CN 

1 


1 kHz 


1.02 kHz ± .13 


56 


-62 


1.5 


2.01 ± .43 


75 


-81 


1.9 


2.24 ± .37 


93 


-100 


2.3 


2.50 ± .39 



The data of Table 2.1 are shown in Figure 2.8. The 
size of the error bars of the experimental data are due to 
the inaccuracies in reading the frequencies of the nodes . 
Note, however, that in all cases the experimental data 
agreed with prerdicted values of the frequency shift with- 
in the accuracies of measurement, so on a qualitative 
level, at least, the theoretical predictions of the nodal 
shifts are supported. 

In further comparison of Figures 2.6 and 2.7, it can 
be seen that although the pattern has been shifted in 
frequency, there is no change in the amplitudes of the 
maxima and only a small change in some of the minima of 
the frequency response pattern when the plate is 
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Figure 2,6, Frequency response, free boundary 
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introduced. This suggests that there is insignificant 
change in the amplitude of the reflected signal# as pre- 
dicted by the reflection coefficient R having a magnitude 
of one. From this# it is inferred that there is insignif- 
icant translation of the incident wave into plate wave 
modes other than longitudinal. The ray-method used in de- 
riving the reflection coefficient R assumes that the com- 
pressional wave incident at the plate boundary remains a 
comressional (or longitudinal) wave in the plate. 

A change in the mode of the wave from compress ional 
to# for instance# flexural or shear# would involve a loss 
of energy in the reflected wave# which would be seen as a 
change in the amplitude of the maxima and minima of the 
frequency response pattern. Since very little change was 
observed# it appears that the incident wave remained a 
compressional wave as it traversed the plate and re- 
entered the fluid below. 

Thus# it has been demonstrated that the effect of 
placing a plate into a fixed-geometry method-of-images 
problem is to introduce only a phase delay into the re- 
flected signal# causing a related upward shift in the fre- 
quencies of maxima and minima in the pressure for the 
fixed source and receiver positions. - 
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Figure 2.8 Frequency Shift due to Plate 



Theoretical 

Actual 



51 



III. RADIATION FROM A POINT DRIVEN 
FLUID-LOADED PLATE 



A. INTRODUCTION 

Chapter II treated the effects on the sound field 
from a source in the water when a plexiglass plate is 
placed on the surface. This chapter will deal with the 
acoustic radiation from a plate driven by a source 
attached to it. When the source is driving the plate, it 
excites the "acoustic mode" of sound propagation in the 
plate, as described in Chapter I. From the previous 
discussion, the acoustic mode leads to radiation away 
from the plate caused by the spatial gradient of the 
amplitude of the flexural wave excited by the source. 

The principal references for the theoretical de- 
scription of the problem are Feit [Ref. 7] and Ross 
[Ref. 21]. A computer program was written to compute 
the radiation pattern from the driven plate at various 
frequencies . An experiment was conducted to obtain 
actual radiation patterns from a plexiglass plate placed 
on the water . 

B. THEORETICAL DEVELOPMENT 

As was seen in the background section of Chapter 
One, the theory concerning the point-driven fluid-loaded 
plate is well known, and its full development will not 
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be given here. However, for continuity, key points 
leading to an understanding of the solution will be 
summarized . 

1 . Radiation from a Propagating Flexural Wave 

A flexural wave propagates in a plate with wave 
number kf where 

kf = -r- (3.1) 

1 °b 

When this plate is in contact with a fluid, vibrations 
of the plate will be imparted to the fluid, producing 
waves in the fluid of wave number k Q where 

k Q » k x + k, (3.2) 

where k x and k ( are the wave numbers perpendicular and 
'parallel to the plane of the plate, respectively. 




Figure 3.1 Sound Radiation from a Flexural 

Wave of a Plate 
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The length of this vector is given by 



V 2 — v 2 .l v 2 

k q *i + K i 



(3.3) 



the component parallel to the plate is 




and from Snell's Law this must be equal to kf • This 
leads to [Ref. 20, p. 163] 




From the boundary conditions 
normal components of particle 



(3.5) 

at the interface, the 
velocity must be equal. 



v lr = u l. 



(3.6) 



where Vj_^ refers to the perpendicular component of the 

plate particle velocity and u L is the fluid particle 

o 

velocity perpendicular to the plate. 

From these conditions, it can be shown that [Ref. 21 
p. 163] the acoustic pressure £ on the plate due to 
the fluid is: 



H “ ” j w Po c o 



V*o 2 - k f 2 



w 



(3.7) 



where w is the normal displacement of the flexural 
wave. Equation 3.7 gives the acoustic pressure £ on 
the plate attributable to fluid loading. Note that £ is 
dependent on the relative values of k Q and k f . 
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This leads to the important effect of coincidence. 
When kf>k Q / the expression for the fluid pressure is 
real, which represents mass-loading [Ref. 21 p. 163] and 
acoustic power is not radiated into the fluid. When 
= kf , the condition of coincidence is reached. The 
expression for the radiated acoustic pressure from the 
flexural wave would seemingly predict infinite pressure. 
(This cannot, however, occur, because fluid loading pre- 
vents the two speeds from being exactly equal [Ref. 21 
p. 163]) When kf<k Q , the expression for the fluid pres- 
sure in Equation 3.7 is imaginary, so the effect is that 
of a radiation resistance, [Ref. 21 p. 163] and the pro- 
pagating wave radiates power away from the plate. 

The coincidence frequency is obtained by equating 
the wave numbers of the flexural wave in the plate and 
the compressional wave in the fluid and solving for fre- 
quency w c . But because, the flexural wave speed varies 
with frequency and formulae for the flexural wave speed 
depend on whether thin plate or thick plate approxima- 
tions are used, there are different forms for the solu- 
tion for the coincidence frequency. 

As there were three approximations for the flexural 
wave speed (Equations 1.6 - 1.8), there are three 
formulae for the coincidence frequency. In the low 
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frequency limit, or thin-plate theory, the coincidence 
frequency is given by: 



o>ci 



12 

t 




(3.8a) 



In the intermediate frequency range, the coincidence 
frequency cannot be easily solved for implicitly, but is 
given by the solution of : 






Cl 



n T 



, “ci .1/2 w ci 

+ — ' 2 - 2— 

«p tbp— 

w. 



, 'Cl. 2 

1 " 




(3.8b) 



In the high frequency range, the coincidence frequency 
can be given by: 




(3.8c) 



Note here that if Cp < Vf c Q , Equation 3.8c is 

negative, which implies that coincidence cannot be 
reached. Recall that the high frequency asymptote of 
the flexural wave speed is the Rayleigh wave speed, a 
value somewhat less than the shear wave speed. If this 
speed is less than the speed of sound in the fluid, the 
condition of k Q = kf, or coincidence, cannot be reached. 
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2. Radiation from a Point Source 



Heretofore, we have discussed the flexural wave 
itself radiating power from the plate. The conclusion 
that no power is radiated below coincidence applies only 
to infinite plane bending waves on an infinite plate. 
If excited by a localized force even below coincidence, 
sound will be radiated because the waves in the plate 
have divergence [Ref. 20 p. 169]. 

The pressure at a field point ( R, 9) from a 
point-driven fluid-loaded plate can be, in a derivation 
by Feit [Ref. 7, p. 1492], eventually expressed in the 
form (3.9) 



£(R, <}), 9 ) = 
where 




(k 0 Rsin<j>sinw) e^o Rs ^ n< ^ cosw V( w) sinwdw 



p(R,<j>,9) = acoustic pressure in spherical 

coordinates with the source at the 
origin. 



R = the slant distance from the source, 

<j> = the azimuthal angle, 

9 = the circumferential angle measured in 
the plane of the plate, 

w = the variable of integration, 

k Q - the wave number in the fluid, 

C = contour of integration shown in 
Figure 3.2, 

F 0 — the amplitude of the source. 
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H = the zero-order Hankel function of the 

° first kind# and 

V(w) = a complicated formula that relates the 
vibrating field of the plate to the water 
through the Mindlin-Timoshenko thick 
plate equation. V(w) contains the poles 
of the integral in Equation 3.10. 

The solution of the integral was obtained by 
Brekhovskikh [Ref. 18, pp. 242-258] and is given by Feit 
[Ref. 7, Equation 24]. Ross, however, gives a more use- 
ful form of Feit ' s solution: 




where Tp is the plate shear parameter, given by 

r p - 2.65(1+ | o + | a 2 ) (3.11) 

and 0 is the fluid loading factor, given by Equation 
3.13, which is discussed fully in the next section. 

The bracketed term in Equation 3.10 along with 
the cos<$> can be described as a directivity function 
H(o>,<t>). This directivity function was encoded in 
FORTRAN and the program (given in Appendix A) was used 
to generate plots of theoretical directivity patterns 
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at several different frequencies. From these figures 
(Figures 3.5-3.11) the coincidence effect can be readily 
observed. The figures were based on thin plate theory, 
and cover frequencies from below to above coincidence. 
They represent a hypothetical plate and are presented 
for a discussion on the effects of coincidence and fluid 
loading. They will be discussed in section El of this 
chapter . 




Figure 3.2 Integration Path of Integral of Equation 3.9 

(from Feit) 

C. THE EFFECT OF FLUID LOADING 

Although Equation 3.10 is valid for all ranges of 
plate and fluid parameters and frequencies, the expres- 
sion may be simplified greatly and further insight may 
be obtained by considering the effect of fluid loading. 
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The fluid loading factor, 0 , 
specific acoustic impedence of the 
reactance per unit are of the plate 
166] defines it as 

„ _ Po c o 

0 “ +■ 

(Dpp t 



is the ratio of the 
fluid to the mass 
. Ross [Ref. 20 p. 

(3.12) 



where 



p Q * density of the fluid 
c Q = speed of sound in the fluid 
a) = radian frequency of the flexural wave 
Pp = density of the plate 
t = thickness of the plate 
Applying the formula for the thin plate coincidence 
frequency, Equation 3.8a, the fluid loading factor 
becomes 



1 Pq c p ^c 
Vl2 Pp c o “ 



(3.13) 



where 

Cp = the plate longitudinal wave speed, given in 
Figure 1.1a, 

uj c = the coincidence frequency. 

A fluid loading factor much less than 1 corresponds 
to light fluid loading, which occurs, for example, in a 
metal plate in air. A fluid loading, factor much greater 
than 1 corresponds to heavy fluid loading found, for 
example, a plexiglass plate on water. The effect of 
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fluid loading can best be seen by considering four 
general cases: frequency below coincidence/ light and 
heavy fluid loading; frequency above coincidence/ light 
and heavy loading . 

When oj<oj c and 0<<1, a case of light fluid loading, 
the radiation pattern is that of a small baffled piston, 
or monopole-type radiation, non-directive. For this 
case. Equation 3.10 reduces to [Ref. 7 p. 1492] 

p(r, <j>, 0 ) = -° - i r e^ koR (3.14) 

2 irppt 

Ross [Ref. 20 p. 170] states that the radius of this 
piston is the sum of the radius of the surface over 
which the force is applied and a quarter of the flexural 
wavelength. Ver and Holmer [Ref. 12 p. 278] give a 
slightly different form. They state that for a plate 
driven by point source, the equivalent radius of the 
baffled piston is 0.286 X c • where X c is the flexural 
wavelength at coincidence. 

When, o}< oj c and 0>>1 , heavy fluid loading, the point 
force acting on a thin plate radiates sound as through 
the plate were not there, representing dipole- type 
radiation, and the radiation pattern is cosinusoidal. 
Equation 3.10 reduces to [Ref. 7 p. 1492]: 

p(R/4>) = ° e j koR cos<|> (3.15) 
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When o)>u) c and £<<1 , the radiation pattern is still co- 
sinusoidal [Ref. 21 p. 171]. 

When a)>a) c and 0>>1 , Equation 3.10 can be reduced 
only slightly, and gives the interesting radiation pat- 
terns corresponding to the flexural wave radiation and 
radiation from the driving piston. 

D. THE EFFECT OF STRUCTURAL ABSORPTION 4 

Above coincidence, then, the absorption of the flex- 
ural wave as it propagates through the plate will effect 
the radiation pattern. Feit [Ref. 7, p. 1494] states 
this and states that (without proof) the far field pres- 
sure maximum is reduced by 5 dB for a loss factor, 
n=0.1. (See Equation 3.18 for the definition of the loss 
factor n .) Ungar [Ref. 22, p. 453] lists a loss factor 
typical of plexiglass of n = 0.02 to 0.04 in the audio 
frequency range. Therefore, it is expected that the de- 
crease in predicted peak pressure amplitude due to the 
absorption of the flexural wave will be less than 5 dB. 

An experiment was conducted to quantify the loss 
factor for the plate used in the experiment. This was 
done also to investigate the effects of the reflection 
of the flexural wave from the edge of the plate, i.e.. 



^Although the effect is the same, the term absorption 
is perhaps more correct than damping as the term damping 
implies additional materials applied to the plate. 



_ 62 



-dieter Va<s ‘ uc e p, 

tUde « ^i ou * d USed *° -a SUre th Piate -P y o, 
** -o US • dista — r avay fr ; he *■* in ^ 

Was d «e to th SUItS Wer e dif fl 

° mster - ° f data **. 

Piats - However, Vari °u* test '. ”** acce ier- 

acc eler 0 jneter Was ^ °° UWin 5 betwe en 31009 ^ 

COUWinP between - obt^ ^ 

<3e stru P f • be data ng the 

t-ive tesf- ^surem Sa ^e 

st va s urej nent s ,• 

A Pr 0per p ° s sibl e Usi a non- 

“ *“■ ■»<- ?r ,h ^ :r 

holt the ° the plat^ v °uid have k 

. the acce -lerometer , at «»• test . ^ 

egu aI tor quin ° 4116 t. at P ° ints ^nd 

Were desir e( j, s 9 ' H ° Wev er, only * tSSt Points 

ai * e P° Ss ibie. ° tMS " as ^t done '••uxt. 

' bln ’ 

°Pti cai _ or smalie 0 U meas urenie n< . 

ac °ustic m eth r s ^Ples of ents ' 

An att ho <?s. ^teriai a , 

at tempt VaQ ' a nd 

a vas ^ 

s Jnaj..£ n • 7^ 

Siraiia P to 7h ° f PleXi9iass < Of S ffl Chiadni Patte rns on 

Saj npie a in the pr operf-v 

^PPr° xijnate ex Periment S 

Iou ^speav y 10 CJn by o n * 'The 

pe *fcer drivo y 2 ° Cfll ) va<* al1 

at fre n Veii in SUs Pended aK 

equ enc ies in ex ce Ss of d ab °v e 

nioda? ^ a bo v e ^ , ■ L ^° dB 

o, w ; m - 

Wo Sohiadni 



63 



patterns were induced. (Other samples, for instance 
of metalj did produce Schladni patterns using this 
technique . ) 

Qualitatively, the effect of the plexiglass loss 
factor n was very apparent. Refer to Figure 3.3 which 
is a plot of flexural wave attenuation, after correction 
for cylindrical spreading. As the accelerometer was 
moved farther away from the source, the corrected ampli- 
tude of the flexural wave steadily decreased with range, 
dropping 40 dB by the end of the plate for the frequen- 
cies of interest. No reflections of the pulse off the 
edges were observed. 

Ver and Homer [Ref. 12, p. 276] give the losses in 
power radiated from the plate due to absorption to be 
e -nkf R/2 

where kf = flexural wave number (3.18) 

R = range from rource 

n = loss factor of the plate material 

using this, with the 100 kHz absorption curve (the only 
one straight enough to be trusted) yields n ~ 0.03 
which is consistent with the values of Ungar [Ref. 21, 
p.. 453]. 

According to Junger and Feit [Ref. 13, p. 249], the 
effects of structural damping on the far-field pressure 
pattern are accounted for, for small loss factors, by 
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Figure 3.3. Flexural wave excess attenuation vs range 



multiplying Cp in Equation 3.10 by (l-jn) l/2 
- (1-1/2 jn)* or u> c by (1+ 1/2 jn). They further state 
that damping plays no significant effect on the maximum 
pressure of the radiation field below coincidence. This 
is also intuitively true because below coincidence the 
propagating flexural wave does not radiate and therefore 
its absorption has no effect. Above coincidence, they 
state (as Feit did in his earlier work) that the effect 
is to reduce the far- field pressure maximum. Figure 
3.4, taken from Feit [Ref. 7, p. 1494], show the effect 
of structural damping on the far-field of a point- 
excited thin plate. 




Figure 3.4, The Effect of Damping on the 

Far-Field Pressure (from Feit) 
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E. 



DISCUSSION OF THEORETICAL PLOTS 



1 . The General Effects of Increasing Frequency 

Figures 3.5 to 3.11 were generated from the pro- 
gram given in Appendix A. They show the radiation 
directivity pattern H ( <f> ) / normalized to the on-axis 
pressure, for a frequency range from below coincidence 
to above coincidence . 

The plate parameters used in these figures were 
those of a hypothetical plate and fluid combination. 
They do# however, show quite well the effects of in- 
creasing frequency on the radiation pattern in the 
fluid. Thin-plate theory was assumed, and the ratio of 
o)/ci) c was the parameter that was varied. 

Figures 3.5 to 3.7 show the directivity patterns 
below coincidence. Far below coincidence , Figure 3.5, the 
radiation from the plate is cosinusoidal, with maximum 
radiation at $ = 0° and minimum radiation at $ = 90°. 
It was shown in Section B1 of this chapter that below 
coincidence an infinite plane propagating flexural wave 
will not radiate sound. It was shown in Section C that 
the radiation pattern should be that of a low frequency 
small radius piston mounted on an infinite baffle. This 
is revealed in these first three figures. Note that as 
the coincidence frequency is approached, the pattern 
changes into a wider (Figures 3.7 and 3.8) pattern. 
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As the frequency of the source becomes larger 
than the coincidence frequency, the phase speed of the 

9 

propagating flexural wave exceeds that of the fluid, and 
therefore the propagating flexural wave can radiate 
energy away from the plate. This is seen in Figures 3.9 
and 3.10. 

Well above coincidence. Figure 3.11, a third 
lobe is predicted. This, according to Feit, is the 
" . . . nonpropagating near-field deflection of the 
plate .... itself becoming a propagating wave." 
[Ref. 7, p. 1493] 

Feit, in a later work with Junger [Ref. 13, p. 
248], state that thick plate theory yields larger values 
for the angles at which the propagating flexural wave 
radiates and smaller values of peak pressure. That is, 
the lobes in the pattern will be at a greater angle and 
have a smaller value. Figure 3.11, from Feit and 
Junger, show a comparison of the size of the pressure 
lobe and the angle at which it appears for thin-plate 
theory and thick-plate theory. The figure neglects the 
effect of structural damping. 
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Figure 3.6. Theoretical response, a)/a) c = 0.52 
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Figure 3.7. Theoretical responce, o)/g) c = 0.79 
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Figure 3.8. Theoretical response, co/co 
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Figure 3.9. Theoretical response, co/co c = 0.94 
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Figure 3.10. Theoretical response, (o/oj 
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Figure 3.11. Theoretical response, aj/o) c = 1.57 




2 . The Effect of the Fluid Loading Factor/ B 

The plots of the previous section showed the 
general trend in the radiation pattern when frequency is 
increased for a given set of thin-plate parameters. 
However, there are subtle changes that take place in the 
patterns when different parameters are used. The den- 
sity of the plate can be measured with reasonable 
accuracy, and Poisson's ratio should be somewhat similar 
among plexiglass samples. However, Young's modulus and 
the exact effect of Poisson's ratio in the wave speed 
formula may vary. Additionally, Equation 3.10 lends 
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itself better to changing certain parameters than 
others . The parameter chosen for variation was the 
fluid loading factor, 3 , for fixed ratios of cd c / . 
Figures 3.13 to 3.21 are the results of these runs. The 
effect of varying the fluid loading factor 0 as follows. 
Increasing the fluid loading smooths the radiation 
pattern, reducing the peaks and making the coincidence 
effect less dramatic. 

Below coincidence, fluid loading should not have 
a great effect on the radiation pattern. This is seen 
in Figures 3.13 thru 3.15. A change in the fluid load- 
ing factor 0 for a frequency ratio of w c / w = 2 below co- 
incidence has a minimal effect on the radiation 
pattern . 

At coincidence, the effects of fluid loading on 
the radiation pattern can be seen: for light fluid 
loading. Figure 3.16 shows a well-developed radiation 
pattern, but as 3 is increased (Figures 3.17 and 3.18) 
the radiation pattern is less pronounced and the off- 
axis pressure gain not as great. 

Similarly, above coincidence, (u) c /o>) — 0.86, as 
seen in Figures 3.19 through 3.21, the increase in fluid 
loading causes a less dramatic radiation pattern. 

The importance of 3 is noticeable. Even though 
in the computer runs 3 was never "much less than one" 
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Figure 3.13 Theoretical response, 8 = 0.65, (d c /o) 
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Figure 3.14 Theoretical response, /9 = 0.68, (o c /a> 
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Figure 3.16 Theoretical response, 3 = 0.32* 
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Figure 3.17 Theoretical response, 3 = 0.34, co c /^ 
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Figure 3.18 Theoretical response 
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Figure 3.19 Theoretical response 
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Figure 3.20 Theoretical response, 3 = 0.29 
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which was described as the condition for light fluid 
loading, small changes in 3 caused large changes in the 
radiation patterns . 

F. DESCRIPTION OF EXPERIMENT 

An experiment was conducted to measure the radiation 
pattern from a point-driven sheet of plexiglass. See 
Figure 3.22. A l/8th-inch thick sheet of plexiglass of 
dimensions 87cm x 152cm was floated on the surface of a 
freshwater anechoic tank. A small silicon bead was 
applied to the edges of the plate to keep it afloat. In 
the center of the plate a Tonpilz longitudinal vibrator 
of the author's design 5 and manufacture was bolted to 
the plate. The vibrator was driven by an HP-3314A 
function generator operating in the pulse mode. A maxi- 
mum driving output voltage of 30V was applied. 

A schematic cross-section of the Tonpilz vibrator is 
given in Appendix B. Its tail mass was made out of 
stainless steel and its head mass out of aluminum. This 
was to impart as large a signal as possible to the plex- 
iglass sheet while keeping the vibrator's natural 
resonance frequency fairly high. An admittance circle 
of susceptance vs. conductance and plots of admittance 



5 For an excellent description of the design of 
Tonpilz vibrators, see Wilson [Ref. 23, Ch. 6]. 
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and suscept^nce vs. frequency for the vibrator are also 
given in Appendix B. They show the motional resonance 
frequency of the vibrator to be 92.4 kHz# and a Q of 
about 112. 

It is stated that although the motional resonance 
frequency of the vibrator was within the frequency range 
of interest for this experiment, this was not of 
key importance because a swept-frequency mode was not 
used. However, what is of interest in the pulse mode is 
the Q of the vibrator . The Q of the vibrator showed its 
effect in the transient response of the vibrator in the 
pulse mode. The transient response decay time grew when 
the driving frequency approached the motional resonance 
frequency of the vibrator. However, in every data set 
measurement, the pulse length was adjusted to give maxi- 
mum steady state response time for the given pulse 
interval . 

The head mass of the vibrator was about 1 cm in 
radius and was attached to the plate with a thin film of 
silicon grease to ensure positive contact between the 
head mass and plate. The clamp rod of the vibrator was 
used to compress the vibrator against the plate. 

The radiation pattern from the driven plate was mea- 
sured with a Celesco LC-10 hydrophone whose output was 
sent through a Tektronics AM-502 low noise preamplifier 
to a Kikusui oscilloscope. The LC-10 was mounted on an 
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HP 331 4 A 
Function Generator 



Kikusui COS 5060 
Oscilloscope 




Figure 3,22, Driven Plate Experimental set-up 
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angular jig that allowed it to be swept from <j> = -90° to 
90° at constant range from the vibrator. An Endevco 
K920 accelerometer was used at times to probe the plate 
to obtain attenuation and edge-reflection data. 

The pulse mode was used because of a problem of 
electrical pickup between the leads to the vibrator and 
the LC-10 . The electrical pickup could not be removed 
even with the use of grounded coaxial cable and shielded 
connectors . 

The length of the pulse varied with frequency. The 
limiting factors for the pulse length were the pulse in- 
terval and the transient response time of the vibrator. 
The pulse interval was limited because the LC-10 was 
located at a fixed distance from the source and recep- 
tion of a pulse by the hydrophone had to be ended before 
the next began to prevent overlapping signals. The 
transient response time of the vibrator also limited the 
length of pulse used. The received pulse showed the 
transient and stead-state response of the vibrator. 
However, amplitude data could only be taken on the 
steady-state portion of the received signal. This 
further limited the time available in the pulse with 
which to drive the plate in a steady-state condition. 

Both the pulse interval and pulse length were ad- 
justed for each data set measurement to provide the 
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maximum time available with which to drive the plate in 
a steady-state condition. 

The properties of the plexiglass used were deter- 
mined experimentally by Johnson and Denny [Ref. 23.] 
Table 3.1 lists the wave speeds for the three wave types 
they measured. Note that both the shear wave speed and 



TABLE 3 . 1 


Experimental Wave Speeds 


for 


Plexiglass 


C P = 


2240 m/ s 


c s = 


1373 m/s 


c b = 


700-1150 m/s 



therefore the high-frequency asymptote of the flexural 
wave speed are less than the speed of sound of the 
water. Therefore# coincidence cannot be reached. 

G. EXPERIMENTAL RESULTS 

1 . Discussion of Plots 

Several plots were obtained at frequencies from 
50 kHz to 245 kHz with varying degrees of success. As 
was stated earlier# the high frequency asymptote, of the 
flexural wave speed was less than the speed of sound in 
the fluid# so the coincidence effect could not be 
obtained. Still there was a change in the radiation 
pattern from the plate as frequency was increased. 
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Figure 3.23 shows the experimental radiation 
pattern when the plase was driven at f=50 kHz. There 
are several aberrations in the pattern when the deflec- 
tion angle increased above 40° . Two plots were made, 
one with the hydrophone at a range of 0.5 m and the 
other at a range of 1 m . This was done to see if the 
distance from the source to the vibrator had an effect 
on the pattern. It can be seen that it did not, as the 
two curves appear to have quite similar behavior. 
Still, the general cosinusoidal pattern of a point- 
driven plate will below coincidence is seen. 

In Figure 3.24, the aberrations of the previous 
pattern are more pronounced. The. next sets of data. 
Figures 3.25 and 3.26 at 92.3 and 93.5 kHz, respectively 
show the development of a side lobe at an angle of about 
45° . Figure 3.26, where data were taken at angles of 
-90° to 90°, show an asymmetry in the pattern. This was 
seen in a few plots and is thought to have been caused 
by asymmetrical excitation of the plate. This is dis- 
cussed further in section G.3. 

The next three plots. Figures 3.27 to 3.29, at 
frequencies of 100 kHz, 120 kHz, and 150 kHz show the 
development of lobes at about 40°. ' There are multiple 
aberrations in the pattern which may or may not have 
been additional minor lobes. Note that only the pattern 
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Figure 3.23. Experimental response, f = 50 kHz 
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Figure 3.24.. Experimental response, f = 68.9 kHz 



92.3-KHZ 
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Figure 3*25. Experimental response, f = 92.3 kHz 



93.5-KMZ 
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Figure J3.26. Experimental response, f = 93.5 kHz 



100-KHZ 
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Figure 3.27. Experimental response, f = 100 kHz 



120-KHZ 



o 




98 



Figure 3,28, Experimental response, f = 120 kHz 



150-KHZ 
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Figure 3.29, Experimental response, f = 150 kHz 
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Figure 3.30. Experimental response, f = 170 kHz 



194-KHZ 
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Figure 3.31. Experimental response, f = 194 kHz 



196-KHZ 
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Figure 3.32. Experimental response, f = 196 kHz 



225-KHZ 
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Figure 3.33. Experimental response, f = 225 kHz 
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Figure 3.35, Experimental response, f = 246 kHz 



at f = 150 kHz show a clean major and minor lobe 
pattern . 

In the following three sets of plots# figures 
3.30 to 3.32 at f = 170 kHz# 194 kHz# and 196 kHz# the 
well-developed pattern of Figure 3.29 is lost and the 
pattern became more complex# with multiple aberrations 
which may have represented small minor lobes in the pat- 
tern or interference resulting from incomplete and asym- 
metrical excitation of the plate. 

The final three data sets# figures 3.33 to 3.35 
at frequencies of 225 kHz# 227 kHz# and 246 kHz, respec- 
tively# show patterns that are very complex# and as evi- 
dence in figures 3.33 and 3.35# very asymmetrical. The 
side lobes of the pattern are just as large as the on- 
axis lobe. 

2 . Error Analysis 

The method of data collection was as follows. 
For a given frequency# with the hydrophone on-axis (0°) 
the output voltage of the preamplifier was adjusted to 
7.0V. All readings thereafter were immediately corrected 
to dB re 7.0V. The signal-to- noise ratio (SNR) varied 
between data runs because of the output response of the 
vibrator. On the average# however# the SNR was esti- 
mated to be 30 dB at zero deflection angle. 
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Since the data were normalized to the on-axis 
reading, the error increased with decreasing sound pres- 
sure (and hence voltage) levels. Table 3.2 shows the 
resulting dB error for various levels. 



EXPERIMENTAL ERROR 


TABLE 3 . 2 


SPL 


ERROR 




OdB 


re 7V ± 


0.12 


dB 


-5 


± 


0.56 


dB 


-10 


± 


0.84 


dB 


-15 


± 


1.27 


dB 



Angular measurements were quite precise with the 
jig that was used. Angles were precise to ±0.1 degree. 
As precise as the angular measurements were, they were 
not as accurate. While every attempt was made to accu- 
rately position the hydrophone directly underneath the 
vibrator at the start of each data collection period, 
parallax and the water refractive effect hampered accu- 
rate positioning. It is estimated that the angular 
accuracy was ±5 degrees. Techniques to improve the 
accuracy were abandoned because it was reasoned that the 
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the result of any inaccuracies would be an angular 
skewing of a pattern, not a change in the pattern 
itself . 

It would seem therefore, that inaccuracies or 
errors in data collection were not the cause of such 
large aberrations are asymmetries in the experimental 
plots . 

3 . Discussion of Results 

Since the high frequency asymptote of the flex- 
ural wave speed was less than that of the speed of sound 
in the fluid, coincidence was not possible. Therefore, 
the large off-axis pressure maxima associated with the 
driven plates above coincidence were not obtained. This 
may, however, be proper modelling of the ice. Stein 
[Ref. 2, p. 288-289] found in his Arctic ice experiments 
that " . . . at no point is the first antisymmetric mode 

[i.e., the flexural wave] supersonic with respect to the 
water; there is no coincidence frequency above which 
phase speeds are faster than the water wave speed. This 
implies that in the plane-wave limit the first antisym- 
metric mode never radiates and never loses energy to the 
water . " 

After a careful search of recent literature on 
driven plates below coincidence, no theories were found 
that would suggest a radiation pattern other than that 
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of a cosinusoidal pattern. Nor were any theories found 
on the radiation patterns from such acoustically slow 
plates that coincidence was not possible. Therefore, 
theoretical plots, such as those developed for earlier 
sections, were not available. However, some possible 
explanations are offered about the complexities and 
asymmetries of the patterns. 

Since this was a finite-sized plate it could be 
argued that plate modes of vibration were set up, which 
caused nodes in the plate that would affect the pressure 
field in the fluid. However, the pulse mode of opera- 
tion was used, which reduces the likelihood that this 
was responsible for discrepancies. Also, it was found 
(and stated earlier in this chapter) that no reflections 
were observed from the plate edge. Therefore, excita- 
tion of normal modes of vibration in the plate are 
thought to be unimportant for the pulse experiment. 

Another possible explanation is that the hydro- 
phone was, at 0.5 m, in the "near-field" of the plate. 
Feit [Ref. 7, p. 1482] used the asymptotic approximation 
k 0 R>>l in the evaluation of the integral. Equation 3.9. 
At the lowest frequency at which data were taken, 50 
kHz, k Q R equals 106, so it would appear that the far- 
field assumption is valid. 

The complexities of the patterns might also be 
explained by the pulse not being long enough. The pulse 
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interval was limited to 0.30ms, based on the speed of 
the pulse through the water to the receiver and the 
transient response time of the vibrator. Table 3.3 
gives values of travel distance that a flexural wave 
packet would achieve in the plate in 0.30 ms for the 
frequencies used. 



TABLE 3 . 3 


f (kHz) 


c^(m/s) 


d(m) 


50 


725 


0.22 • 


100 


1026 


0.31 


200 


1150 


0.35 



The cosinusoidal patterns expected below coinci- 
dence are based on an assumption of cylindrical diver- 
gence of the flexural wave on an infinite plate. That 
is, the pressure at a field point in the fluid is made up 
of the sum of pressures caused by the deflection of each 
incremental area of the plate. If the entire plate were 
not excited, this sum would be incomplete. From Table 
3.2, a pulse of flexural wave of frequency 50 kHz will 
only travel out 0.22m in 0.3 m/s. Thus, the entire 
plate was not being excited. Even at 200 kHz, . the 
entire plate was not excited. This certainly could be 
one of the contributing factors in the complexity of 
experimental plots. 



Prior to the discovery of the electrical pickup 
problem mentioned earlier in this chapter, the plate was 
driven in the CW mode. The complexities in the patterns 
were still there. This would indicate that use of the 
pulse mode was only a contributing factor, not the main 
cause of the complexities of the experimental plots. 

A more important factor is thought to be the ab- 
sorption of the flexural wave in the plexiglass. 
Following the same argument above, that is, the pressure 
at a field point is the sum of the pressures caused by 
the individual deflections of each incremental area of 
the plate, the effect of the flexural wave absorption 
would be to greatly reduce the contributions of radia- 
tions from the incremental areas far from the source . 

The so-called "acoustic mode", of sound propa- 
gation, which this experiment was modelling, assumes 
sound radiation away from the source in the fluid from a 
lossless plate. Plexiglass is a highly absorptive 
medium; therefore, the propagating flexural wave decayed 
in amplitude due to losses as well as divergence as it 
travelled away from the vibrator. Consequently, the 
incremental areas of the plate far from the source did 
not deflect nearly as much as those close to the course, 
so they did not contribute as much to the pressure at a 
given point in the fluid. 
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The plate might be thought of as a complex aper- 
ture, or a quasi-piston, with a radius of the head mass 
radius and a portion of the adjacent plate. A rough 
estimate of the aperture radius is the head mass 
radius plus either one- fourth the flexural wave radius 
plus either one-fourth the flexural wave radius or one 
plate thickness. Both give similar results. Based on 
this, k Q a for radiation into the fluid varies from 2.79 
for a frequency of 50 kHz to 14 at f = 250 kHz. 

Even though the aperture is shaded and phased, 
it is noteworthy that for a k Q a in excess of 4, a piston 
is beginning to develop a side lobe, and for a k Q a of 
10, it has about two side lobes. This is not inconsis- 
tent with the angular scale of the lobelike structure 
observed . 

The asymmetry seen in a few of the radiation 
patterns might have been caused by radially non-uniform 
excitation of the plate. As was mentioned in section F 
of this chapter, the vibrator was compressed against the 
plate by a bolt. However, asymmetrical compression of 
the vibrator against the plate was possible. Possible 
causes of asymmetrical compression are improper machin- 
ing of the head mass, imperfect centering of the bolt 
through the vibrator assembly, or a non-uniform layer of 
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grease between the head mass and plate. A non-uniform 
radial excitation of the plate seems to be a plausible 
explanation for asymmetry in the radiation patterns . 

Two important points that can be made are that 
it appears that the radiation pattern is a function of 
the source, not any flexural wave propagation. 
Secondly, situations where coincidence is not possible, 
such as a plexiglass plate on water, present radiation 
patterns at lower frequencies similar to radiation pat- 
terns below coincidence in situations where coincidence 
is possible. 
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IV. CONCLUSIONS 



This study was a low-cost, "low-tech", first-cut 
attempt to duplicate a sound source mounted on the ice 
in the Arctic. The first experiment performed showed 
that when a plate is introduced at the boundary of a 
Lloyd's-mirror-type problem, the effect of the plate was 
only to introduce phase delay, causing a shift in the 
node and antinode pattern. If any flexural waves were 
excited in the plate, which might have caused an in- 
crease in amplitude in the received signal, they would 
have been absorbed by the plexiglass. 

In the second experiment performed, it was found 
that even though the theory behind a point-driven fluid- 
loaded plate is well known and many interesting results 
can be gleaned from that theory, extracting experimental 
plots upon which a theory can be based is difficult. 

Plexiglass was originally chosen because its speci- 
fic acoustic impedance is fairly close to both ice and 
water. However, the flexural wave absorption suffered 
in commonly available plexiglass and more importantly, 
its low flexural wave speed render it a poor choice to 
use in the lab in this type of experiment, where wave 
generation and propagation are important. 
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Several interesting types of plastics are 
available. Often the manufacturer lists the density and 
Young's modulus of each in its catalog. However, flex- 
ural wave attenuation data are never available. It is 
recommended that the next step in research along these 
lines be to collect samples of each type of plastic and 
measure these parameters. The time available for this 
study prevented such an analysis. 
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APPENDIX A 



COMPUTER PROGRAM 
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*900010 

*900020 

*900010 



THIS PROGRAM I79LEWXTS ROSS FORMUJ1 4.41. >9104 IS A 
SOLUTION ran TV* PAR-FIELD WKJlil PROM A POINT-DRIVEN 
*.UID-LQACED PLATE AS OR 10 1 HALLY DERIVED IV REIT. 



IT IS A HON- INTERACTIVE PROORAM. 
PARAMETERS. AMD MATO SOUK) VELOCITY 
TV* PROGRAM HUM. 



OR INTEREST. PLATE 
TO M INSERTED INTO 



T>* POOOUN USES 0IS3RLA VERSION 9.0 (JOFTNARS PRODUCT OF IS3C0. 
SAN OICOO. CA) TO PRODUCE A GRAPHICAL OUTPUT. TV* QRARN IS A 
POLAR PLOT WITH TV* I -PI TO PI) AXIS MINS TV* PLATE AIQ A POINT 
SOURCE AT T)« ORIGIN. 

TV* PROGRAM IS CURRENTLY SET PQR A OK-f IONTN INCH 9WZT OF PLEXI- 
B LA ST IN PRCSM4ATER 

TV* PARAMETERS AREi 

X IS 1>* HAVE NUMBER IN TV* FLUID 

PSI. IPSI ARE T>« REAL AHO INTEGER VALUES OF T>* ANGLE 
RPSI IS TV* ANGLE EXPRESSED IN RADIANS 
PI • 1.141ft 

NORM IS THI NORMALIZATION FACTOR. 

P IS TV* FREQUENCY IN HERTZ 
M. OCGA. IS THE RADIAL FREQUENCY 
RMO IS THE DENSITY OF T>« PLATE 
RMONT IS THE DENSITY OF TV* HATER 
M IS T>* TNIOOCSS OF THE PLATE 
N IS T* AREA DENSITY OF T)« PLATE • I 
C IS THI SPEED OF SOUND IN T>* M.UI0 
Cl IS T>* PLATE HAVE VELOCITY 



TV* ON-AXIS PRESSURE 



• R 



Cl - IV/RHOI I — SI®4A»»2 J .1 
PI AM© PI AR* PARAMETERS USED IN TV* PROGRAM i 



PI IS TV* RATIO OF TV* SOM 

fi • ic/cn»»t 

FI IS TV* RATIO OF TV* FREQUENCIES. SQUARED! 

Ft • IH/NOUOC ) **t 

S AM4MP IS THE PLATE $>«AR PARAMCTER • ].«• FOR PLEXIGLASS 

• S .29 FOR STEEL 
i J.J9 FOR ICS 

NSUOC IS TV* RADIAL COINCIDOCE PRTOUENCV 

MOURCH IS TV* HIGH FREQUENCY APPROX IMATIGN OF N SURC 
NOUKL IS TV* LON FREQUENCY APPROX I HAT 1 04 OF NSURC 

01 TV90U94 GS ARE INTER9VDIATE QUANTITIES USED IN TV* PROGRAM, 



REAL X.PSl. RPSI .R. PI. NORM. KPM.I04.C3 .OAMWAP 
. PA. Ft .PC.PO.A.O . 0 . MSU1CH . HSURCL 

oixnsion p(itn.psintn.RPsinti) 



INTER PARAMCTER VALUES XREl 



PI/I. l«lft«/. C/ ./.Cl/ 
/ /.RAO/ /.RHONT/ 



*994090 
*900040 
*900070 
*900090 
*900090 
*900100 
*900110 
*900120 
*944110 
*900140 
*940190 
*900140 
*900170 
*940100 
*900190 
*900200 
*900210 
*940220 
*9002X0 
*9 002(0 
*900210 
*90 0240 
*900270 
*900290 
*900290 
*900X00 
*900X10 
*900X20 
*900XX0 
*900X40 
*900Xt0 
*900X40 
*900X70 
*900X90 
*940X90 
*900400 
*900410 
*900420 
*9044X0 
*940440 
*900410 
*944444 
*944470 
*900440 
*900490 
*900100 
*900110 
*900120 
*9001X0 
*900140 
*900110 
*9Q0S40 
*900170 
*900400 
*900194 
*900400 
*900410 
*900420 
*9 004X0 
*900440 
*900410 
*900444 
*904470 
*940444 
*940490 
*900700 
*940710 
*900720 
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ENTER TVC FREQUENCY HERE (REAL NUMBER. IN NZ. I 



P • 






XRE. TX X LIMIT OP TX FOX HAVE SPEED IS COMPUTED 
AMD TX FLATE WAVE « BASED ON THAT SPEED 
VW • 0 . 92 • Cl/ 3. 

KP • M / VPH 

XM TX PRODUCT K • H IS COXUTEDt 
KM • KP • M 

TX RACIAL COINCIDENCE FREQUENCY IS CALCULATED HERE* 

TX LP COINCIDENCE FREQUENCY IS NSUtCL t 

MEUBCL • ( 12. "*0 .51 • (C--2.0) / (H*CI ) 

MSUBCH IS TX X LIMIT FOR TX COINCIDENCE FREQUENCY* 
WSUBCN • Cl • (S.44/H) • (C”2 ./ (Cl "»2. - 2.S ■ C"»2.)) 



HSUBC • M3UBCL 

CRAT IS TX RATIO CO/CP. ( C/Cl I 
CRAT • C/CI 



IPRCQK - WTCP/I000.) 

THIS AREA SETS IF TX PAOE AND DIRECTS TX OUTPUT* 
CALL SXRPA ('XNOSTK* » 'A* .3) 

CALL TEX41C 
CALL FRTPLT (72.41 
CALL KjOMUP (0.75) 

CALL PAOE (11.4.51 
CALL AACA2D(«..4.! 



THIS AREA PRINTS TX ANQLES AND MESSAGES t 
CALL MESSA4 CP • t* . 100. S.25 .4.5) 

CALL INTNQ (IPRCQK* *ASUT* • * AtUT* I 
CALL RSALNO (P.0. 'ABUT* . ' ABUT* I 
CALL XSSA4 ('HOC t* .100.' ABUT* , 'ABUT* I 
CALL XSSAQ ('•*.! ,3.92.4.9$) 

CALL MESSAO CIO* .2*2. 94*4. ) 

CALL MES3AQ (* 10* .2,4.94.4. I 
CALL HESS AO ('20* .2*1.45.5.4) 

CALL HESSA4 ( *20* .2,5. 95. 5 .4 ) 

CALL XSSA4 ('39* .2.0.94.5.35) 

CALL HESSA4 CIO* .2,4.95,5.15) 

CALL HE33A0 (*44 * ,2*0.00.4.45 ) 

CALL MES3A4 ( *44* .2*7. 94.4.40 ) 



CALL COXL* 

CALL NOCXX 

CALL FOjORO (-20.0) 

CALL XMAX ('H(PSS) 

. 4 *. 100 ) 

CALL XAXCND CNOP1RST* ) 

CALL XSSA4 (*DB4* . 140.0. .-4.S0I 
CALL XI NT AX 

CALL POLAR ( I. .10.4. .0.1 

THIS AREA IS TX COMPUTATIONAL FART OP TX PROGRAM t 
H m I 

M > 2. • PI ■ P 
PS • M/M4UBC 
K • M / C 



>• 8MP007SQ 

>• BMP00740 

• BMPOO7S0 
BHP00740 

• 4X00770 
9X0 0 74 0 

• 4X00790 

>• 4X00400 

4X 00410 

>• 4X 004 20 

■ 4X 004SQ 
4X00440 

■ BX004S9 
9X00440 
4X00470 
4X00444 
4X00490 
4X00900 
4X00910 
4X00920 
4X009SO 
4X00940 
9X00950 
8X00440 
9X00470 
9X00940 
9X00990 
4X01000 
9X0 1010 
9X01020 
BXOIOSO 
9X01040 
4X01050 
4X01040 
8X01070 
1X01080 
4X01090 
4X01100 
4X01110 
4X01120 
4X0US0 
4X01140 
4X01150 
4X0 1140 
4X41170 
4X41140 
4X01194 

. 4X01204 
BXOI2I0 
4X01220 
4X01230 
4X01240 
4X01250 
4XOI240 
4X01270 
4X412S0 
4X01290 
4X41 S00 
4X41S1Q 
4X0 1320 
4X01SSO 
4X0IS40 
4X01S50 
4XOIS40 
4X01370 
4X0 1340 
4X01390 
1X01400 
4X01410 
4X0 1420 
1X01430 
1X0144Q 
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•CTa • (C • WOfT I / CM*WO»H ) 

FI m (C /Cl ) 2 

n • CM / MUK ) •• 2 

DO It IF1I ■ -90. 90 

FtKM) • RJUT(IPSI) ■ Ff / in. 

v«i cm) ' mj* mwi'M ) ■ fi / ito. 

01 • (SIM<FS1(M)))»«2-#1 

ox • ($ iM ( wi ( Mn )■♦ 



F(N)*COS(FtlCN))*(I.M(eaS(FlI(N))/irrA>»((l.-FI»ai«C2)/(I.»Fl 
0AmAF*QI ) ) )••< )••(-.!) 

F(N)» 20. • AU»10(P(H)/NORM(M.C.Cl.M3UK.aAmA^,trrAn 
this fofco thc minimi* fhexsiwi to ic -zo-ot (as a >«lf im 

SCA4.IM0 TV® AXIS) 

IF (FIN) .Lf. -20 .0 ) F(N) • (-20.01 
M - MM 

THIS *m POINTS T>* MCSSAGC A*OUT T>« FLUID LOAOlWi 
CALL 9C3SA4 ('Misuse • 1*. 100. 1.2.1. 71) 

CALL A CALM) ( PS *2. * AtUT* • * AtUT* ) 

CALL WSSAO CtCTA • I*. 100. 1.2. 4. 91) 

CALL AEALMO (SCTa. 2. * AtUT* . ' AtUT* 1 

CALL *SSA4 CC0/C P • «*. 100.1.2.4.21) 

CALL AlALMO (OtAT.2* ' AtVT* • *AR/T* 1 
CALL XXSJ+ ( ** • M • I* * 100.1.2*4.41) 

CALL KtALHO IXM.2 . * AOVT* . ’ABUT* ) 

CALL ■JtfC(S..2.S.2.2.I.0.2.0) 

THIS OCA POINTS OUT TV* CUAVfi 

call turn* 

CALL ■jO*v (0.02.0.0) 

CALL MAMS* (IS) 

CALL TMCOrv (.OS) 

C MX CUIW (ATS I .0. 101.0) 

CALL OtSCT (’TtOCCXY*) 

CALL 00100.1) 

L ( 0 ) 



ITOP 

M 



TMl FOLLWIMt FVCTl®* SUSPOOOOAM CALCULATES TV* MOK9MLI 2AT I9< 
PACT* - Ml TM VALUES AA| HOAmLIIQ TO T»« »-**lS WMUtf 
POO T>* DCS I AO FOCOUOCY. XXX. ALL FIKSSlItCS AT FS1 ■ 0 DO 



FVnCTIOi MOM (M.C.C1.NOUSC. OAMUF.OCTA) 
01AL FS1 . FI. MOM 

DATA FI /S.141S9/.FSI /•./ 

FI • (C /Cl ) •• 2 

FI • (M / MOUSC ) •• 2 



01 • (ttM(FSI ))»»t-#I 
01 • ( t iMIPtl ))•• 
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•*01710 
09*01740 
09*01770 
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•*01940 
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•*01940 
•*01970 
■9*01900 
■9*01990 
•*02000 
•*02010 
■9*02020 
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APPENDIX B 



DATA ON TONPILZ VIBRATOR 




Figure B.l. Diagram of Tonpilz vibrator 




A: G B: B o MKR 92 540.000 Hz 

A MAX 2.500 mS REAL 975.586 jjS 

B MAX 1.000 mS IMAG -174.457 jjS 




Figure B.2. Susceptance and Conductance vs. 

Frequency of Tonpilz vibrator 
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o MKR 
REAL 

991.492 mS 
I MAG 

178.144 ixS 



<Hor izontal > 

A: G 
A MAX 
1 . 400 rnS 
A MIN 
0.000 S 

< Vertical > 



B: B 
B MAX 
1 .000 
B/DIV 
140.0 



92 200.000 Hz 



mS 



MS 




TONPILZ VIBRATOR ADMITTANCE CIRCLE. 



Figure B.3. Complex Admittance circle of Tonpilz 

vibrator 
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